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The KTeV E799 experiment has conducted a search for the rare decays Kl — > 
and Kl — > 7r'V V — s> 7r°7r° fi + fx~ , where the X° is a possible new neutral boson that was 
reported by the HyperCP experiment with a mass of (214.3±0.5) MeV/c 2 . We find no evidence 
for either decay. We obtain upper limits of Br(A' L ->■ tt°tv°X -> iy°tt° (i + < 1.0 x KT 10 
and Br(K L -> ttVVV") < 9.2 x 1CT 11 at the 90% confidence level. This result rules out 
the pseudoscalar X° as an explanation of the HyperCP result under the scenario that the dsX° 
coupling is completely real. 

PACS numbers: 13.20.Eb, 13.25.Es 



The HyperCP collaboration has reported the possible 
observation of an X° boson of mass (214.3±0.5) MeV/ c 2 
decaying into ll + li~ based on three observed events in 
a search for the decay S + — > pn + n~ [1]. The con- 
fidence level within the Standard Model for all three 
events to overlap within the HyperCP dimuon mass res- 
olution of 0.5 MeV/c 2 is less than 1%. As the X° would 
presumably be a strange-to-down neutral current, it is 
natural to look for it in the kaon sector, specifically in 
the mode K L -> 7r 7r°V° -> 7r°7rV + M~- This letter 
presents the first attempt to detect the rare decay modes 
K L -> 7r°7rV+/i- and K L -> 7r 7r°V° ->• 7r°7r°/i + /^- 

Using a two-quark flavor changing coupling model in 
which the X° couples to ds and theoretical es- 

timates of the K L —> tt 7t X —5- 7r°7r°/i + /i _ branching 
ratio were determined for a pseudoscalar X° and an ax- 
ial vector X° [2]. Reference [2] uses the known value 
Br(AT ± -)• 7r ± /i+/i-) = 8.1 xKT 8 [3] to rule out the 
possibility of a scalar or vector X° as explanations of 
the HyperCP anomaly. These predictions assume real 
dsX° couplings, gp; for a complex coupling with a dom- 
inant imaginary term, |9?(<?p)| > 0.98|<?p|, the predicted 
upper limit is much smaller [4]. Another prediction of 
Br(Ki — > 7r°7r V° — > 7r°7r°/! + /i~) for a pseudoscalar 



X° has been made [5]. Finally, the branching ratio for 
Kt, — > ir Q ir°X° — > 7r°7r°77 has been estimated using a 
sgoldstino model [6]. These results are summarized in 
Table I. 

The E391a collaboration has reported [7] an upper 
limit Br(K L -)• 7r 7r°W° -> 7r°7r°77) < 2.4 x KT 7 , which 
rules out the sgoldstino model of this decay. The pos- 
sibility [8] that X° could be a light pseudoscalar Higgs 
boson of the next-to-minimal supersymmetric Standard 
Model (NMSSM) was investigated at e + e~ colliders by 
CLEO [9] and BaBar [10,11,12] and at the Tevatron (DO) 
[13]. No evidence for a NMSSM light pseudoscalar Higgs 
boson was found. 

The Kl —5- irirX modes have an extremely lim- 
ited phase space. The phase space of Kl — > Tr -^ ^ 
is ten times larger than the phase space available to 
Kl — > n + ir~X a , motivating the search for the former 
over the latter. We have searched for Kl — > 7r°7r°/x + /x _ 
and K L -> 7r 7r°V° -> 7r°7rV+^- in data from the 1997 
and 1999 runs of KTeV E799 II at Fermi National Accel- 
erator Laboratory. 

The KTeV E799 experiment produced neutral kaons 
via collisions of 800 GeV/c protons with a BeO tar- 
get. The particles created from interactions with the 
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X° -»■ /xV Model 


Br(^L ttVX ) 


Pseudoscalar ($t(gp)) [2] 


(8.3l^)xl0" 9 


Axial Vector (5%a)) [2] 


(l.O^)xlO- 10 


Pseudoscalar (|3(ffp) > 0.98|gp|) [4] 


< 7 x 1CT 11 


Pseudoscalar (K((/p)) [5] 


8.02xl0~ 9 


sgoldstino (X° -> 77) [6] 


1.2xl0~ 4 



TABLE I: Summary of predicted branching ratios for Kl 
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target passed through a series of collimators, absorbers 
and sweeper magnets to produce two nearly parallel Kl 
beams. The Kl beams then entered a 65 m long vacuum 
tank, which was evacuated to 1 ^Torr. A doubling of the 
spill length and an increase in instantaneous luminosity 
between the 1997 and 1999 data-taking periods resulted 
in a factor of 2-3 increase in per-spill protons on target 
in 1999. A 10.3 inch beryllium absorber was introduced 
in the 1999 data-taking period to reduce neutron back- 
grounds. 

Immediately downstream of the vacuum region was a 
spectrometer composed of an analysis magnet between 
two pairs of drift chambers. The momentum kick im- 
parted by the magnetic field was reduced from the 1997 
value of 0.205 GeV/c to 0.150 GeV/c in 1999 to increase 
the acceptance for low momentum charged particles. The 
momentum resolution of the spectrometer in 1997 was 
np/P = 0.38% © 0.016%P [14], and in 1999 the momen- 
tum resolution was a P /P = 0.52% 0.022%P. 

The electromagnetic calorimeter was constructed of 
3100 pure Csl crystal blocks arranged into a 1.9x1.9 
m 2 array. Each Csl crystal was 27 radiation lengths 
long. Two holes were located near the center of the 
calorimeter to allow for passage of the beams. The 
electromagnetic calorimeter had an energy resolution of 
<7 E /E ~ 0.4% 2%/^/E[GeV] and the position reso- 
lution was about 1 mm. The muon ID system used a 
Pb wall, three steel filters and three scintillator counter 
planes to identify muons by filtering out other charged 
particles. The muon ID system contained 31 nuclear in- 
teraction lengths of material and had a charged pion fake 
rate of (1.69 + 0.17P [GcV/c]) x 10~ 3 , where P is the 
track momentum. A photon veto system detected pho- 
tons outside the detector acceptance. The upstream sec- 
tion of the photon veto system had five lead-scintillator 
counter arrays located inside the vacuum decay region. 
The downstream section of the photon vetos had four 
lead-scintillator arrays that framed the outside of the last 
three drift chambers and the Csl calorimeter. A more de- 
tailed description of the KTeV detector and photon veto 
system can be found in [15, 16, 17]. 

The signal modes and normalization mode (Kl — > 
7r^7r^7r^) , where one photon was lost down the beam 
hole and tt° d — > e + e~j) were collected by different trig- 



gers. The triggers required in-time energy clusters in the 
calorimeter of at least 1 GeV. The signal mode required 
one (two) such clusters for the 1997 (1999) data-taking 
periods. Two hits were required in the two most down- 
stream stations of the muon system in 1997; in 1999, the 
number of hits required in the middle station was reduced 
to one. The normalization mode trigger required at least 
four in-time clusters and two tracks. 

Both tracks were required to form a good vertex within 
the vacuum decay region, to match (within 7 cm) a clus- 
ter in the Csl calorimeter and to deposit less than 1 GeV 
of energy in the Csl calorimeter, consistent with a muon 
hypothesis. 99.9% of muons with a track momentum over 
7.0 GeV/c satisfied the last three requirements. Each 
of the three scintillator counting planes in the muon ID 
system were required to register at least one hit. The 
invariant /x + pT mass, M^, was required to be less than 
0.232 GcV/c 2 , which is slightly above the kinematic limit 
given by M K - 2M T . 

Four clusters in the calorimeter without associated 
tracks were required. The resolution of the z-vertex de- 
termined from the two 77 vertices associated with a 7r°7r° 
was better than the resolution of the z-vertex from the 
two muons. We considered each possible 77 pair to find 
the combination with the best agreement between the 
positions of the two 77 decay points under the hypoth- 
esis that each originated from a ir° decay. A minimum 
pairing chi-squared, x 2 , was calculated to determine the 
best agreement between the positions of the two 77 decay 
points. A weighted average of z-vertex values for each 77 
in the pairing with the minimum xl was used as the de- 
cay vertex for the event. This vertex was then required 
to be located within the length of the vacuum decay re- 
gion. A 77 mass, M 77 , was calculated for the event using 
the decay vertex from the minimum \ 2 Z pairing. M 77 was 
required to be within 0.009 GeV/ c 2 of the n° mass. 

The Kl — > 7r°7r°/i + /!~ simulation was modeled as a 
four body decay using a constant matrix element. The 
Kl — > 7r°7r°A' n — »■ 7r°7r simulation was modeled as 
a three body decay with a flat phase space, where the X° 
underwent a prompt decay to 

The signal regions for the 1997 and 1999 data were 
based on the M AtAt7777 , Pri^^llll) an d |Pt(mm) — 
Pr(llll) \ resolutions calculated in the simulation. Here 
p\ is measured transverse to the direction of the Kl, 
determined by the line connecting the BeO target and 
the vertex. For a well-measured decay, p\ (^7777) and 
~ Pt(7777)I should be close to zero. The Kl 
signal region for the decay Kl — > 7r°7r°/Lt + /x _ was de- 
fined as 0.495 GeV/c 2 < M W7777 < 0.501 GcV/c 2 and 
p\ < 1.3xl0~ 4 (GcV/c) 2 . The X° signal region for 
the Kl — > 7r°7r°A' — > 7r°7r°/i + /i~ decay was defined 
as 213.8xl0~ 3 GcV/c 2 < M w < 214.8xl0~ 3 GeV/c 2 
and |j4(mm) -Pt(7777)I < 7.0xl0- 4 (GcV/c) 2 . The 
bound on was determined from the conservative hy- 
pothesis that the observations made by HyperCP reflect 
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FIG. 1: a,b) p% vs. M W7777 for the 1997 and 1999 K L -> 
7r°7r° n + fi~ simulation respectively. The boxed signal region 
contains 90% ol all events. c,d) \p^(/j,fx) — Pt(7777)| vs. Af MA1 
for the 1997 and 1999 K L -)■ 7r°7r X° -» 7r°7rV + M~ simula- 
tion respectively. The boxed signal region contains 95% of all 
events. The four plots are shown after all analysis require- 
ments were applied. 



the natural width of the X° [18]. Figure 1 shows p\ 
vs. invariant mass plots from the Kl — > 7r°7r°/i + /i~ and 
Kl — > 7r 7r°X° — > 7r 7r°/z + /x _ signal mode simulations. 



Every Kl decay mode with two minimum ionizing 
tracks and at least one photon was considered as a po- 
tential source of background. Accidental time-coincident 
activity created from particle interactions in the vacuum 
window, neutrons from the target, cosmic rays, beam 
interactions or another kaon decay in flight can over- 
lap with the primary kaon decay in an event to repro- 
duce the signal mode topology. We have simulated all 
known backgrounds to the extent possible. Accidental 
activity was included in the simulation of all background 
mode events. Small branching ratio backgrounds such 
as Kl — ► tt°tt fM^Vp [19] and Kl — > n + ir~j were simu- 
lated with the full statistics of the data. Large branching 
ratio modes such as Kl —> tt p^v^ were also studied 
extensively, although simulated samples with statistics 
similar to the data were not feasible. We find that when 
accidental activity reproduces the signal mode topology, 



Pt 



and the invariant mass M, 



are pushed to val- 



ues well above the signal region. The conclusion that the 
background is negligible is confirmed in the data. 

The normalization mode shares the topological trait 
of four photons and two tracks with the signal mode and 
has a well-understood branching ratio. The vertex in the 
normalization mode analysis was required to be located 
within the vacuum decay region. The signal region for 
the 1997 data was defined by 0.494 GeV/c 2 < M ee7777 < 
0.501 GeV/c 2 andp|(ee7777) < 0.00015 GeV 2 /c 2 . The 
signal region for the 1999 data was a contour that was 
derived from a joint probability distribution of M ee7777 
and P5"(ee7777) signal resolutions from simulations. 

The flux, Fk, obtained from the normalization mode, 
is the estimated number of Kl decays in the vacuum 
decay region. Uncertainties in Fk originated from the 
branching ratio used to calculate Fk and the muon ID 
system efficiency. The uncertainty due to normaliza- 
tion mode requirements was studied by varying the se- 
lection requirements in simulation and data and noting 
the change in the estimated flux. The uncertainty due 
to differences in the signal and normalization mode sim- 
ulations was estimated by varying the selection require- 
ments of both modes in the simulation. The statisti- 
cal uncertainties on the signal mode simulation were less 
than 0.14% for each decay mode. The statistical uncer- 
tainty for the normalization mode simulation was less 
than 0.37%, while the statistical uncertainty for the nor- 
malization mode data was less than 1.14%. Systematic 
uncertainty in the muon ID efficiency came from model- 
ing of the energy loss in the muon filters and from sim- 
ulation of gaps between scintillator paddles in the muon 
planes [20]. A systematic uncertainty associated with 
the muon trigger inefficiency was determined by selecting 
clean K^ decays from a minimum bias trigger [21]. Re- 
sults from these systematic uncertainty studies are given 
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Systematic Uncertainty on Fk 


AF K 
Fk 


VdlldllUll Ul 1 \ Ul LllCLLlZiCLb LKJ1L 

Requirements 


3.57% 


Variation of Signal/Normalization 
Requirements 


5.35% 


Muon Trigger Inefficiency 


2.00% 


Cracks in Muon Counting Planes 


0.50% 


Energy Loss in Muon Filters 


0.40% 


Br{K L -> ttW) [3] 


0.61% 


Br(TY° -> 77) [3] 


0.03% 


Br(Tv° -)■ e+er-y) [3] 


2.98% 


Total Systematic Uncertainty 


7.41% 



TABLE II: Summary of systematic uncertainties on the ap- 
parent Kl flux, labeled as Fk- 



in Table II. 



Br{K L ir°ir°X a ->• ttV/z+aO < 1.0 x lO" 10 . 

Our result for Br{K L ->■ ir°ir°X ->• ttVV+aO is 
nearly two orders of magnitude smaller than the expected 
branching ratios for Kl — > Tr Q ir°X° — > 7r°7r°/i + /!~ from 
[2] and [5], in which X° was taken to be a pseudoscalar. 
This rules out the pseudoscalar X° as an explanation 
of the HyperCP result under the premise that gp is 
completely real and also places a tight bound on gp of 
I^(<7p)I ~ 0.98|gp| [4]. Finally, our upper limit challenges 
the axial- vector X° explanation of the HyperCP result. 
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FIG. 2: Pt vs. My. 



plot for the combined 1997 and 



1999 data sets, which are indicated by triangles and squares 
respectively. The signal box is open. 



The 1997 (1999) signal mode acceptance was 3.14% 
(4.03%) and 2.80% (3.74%) for K L -> 7T°7rV + M" and 
K L -)• 7r°7r°A: -> 7r°7rV + M~ respectively. The 1997 and 
1999 normalization mode acceptances were 4.21 xlO -6 
and 3.26 xlO -6 respectively [22]. While signal mode 
acceptance would drop in scenarios for which the X° 
does not decay immediately, it would not be sharply 
reduced for values of cr < 3 mm, where r is the 
F K was 2.73xlO n in 1997 and 
The single event sensitivity was 
and 4.34X10" 11 for 
Figure 2 displays the 
results of the blind analysis; no events are inside the sig- 
nal regions after opening the signal boxes and no events 
were found within the available [i + fJ.~ phase space. Us- 
ing the method of [23], the 90% confidence level upper 
limits are Br(K L -> 7r°7rV + A*~) < 9.2 x 10~" and 
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4.12x10" in 1999. 
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